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Abstract: Diffuse emission from the Millcy Way dominates tiie 7-ray sky. About 80% of tlie higli-energy 
luminosity of the Milky Way comes from processes in the interstellar medium. The Galactic diffuse 
emission traces interactions of energetic particles, primarily protons and electrons, with the interstellar 
gas and radiation field, thus delivering information about cosmic-ray spectra and interstellar mass in 
distant locations. Additionally, the Galactic diffuse emission is the celestial foreground for the study of 
7-ray point sources and the extragalactic diffuse 7-ray emission. We report on the latest developments in 
the modelling of the Galactic diffuse emission, which will be used for the Gamma Ray Large Area Space 
Telescope (GLAST) investigations. 
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Introduction 

Diffuse Galactic emission (DGE) dominates the 7- 
ray sky with more than 80% of the total luminosity 
coming from processes in the interstellar medium 
(ISM). The DGE is a tracer of energetic interac- 
tions of cosmic ray (CR) particles in the ISM, and 
is produced by inverse Compton scattering (IC), 
bremsstrahlung, and 7r°-decay. It delivers infor- 
mation about spectra and intensities of CR species 
at distant locations and allows the study of CR 
acceleration in sources as well as propagation in 
the ISM. Gamma rays can be used to trace the 
interstellar gas independently of other astronomi- 
cal methods, e.g., the relation of the molecular H2 
gas to CO [19] and hydrogen overlooked by other 
methods [7]. Additionally, the DGE is the bright 
"background" against which 7-ray point sources 
are detected and its accurate determination is im- 
portant for localisation of such sources and their 
spectra, especially at low Galactic latitudes. Fur- 
thermore, the DGE acts as a "foreground" for any 
extragalactic signal that we seek to recover 



Calculation of the DGE requires a model of CR 
propagation. Such models are based on the theory 
of particle transport and interactions in the ISM as 
well as many kinds of data provided by different 
experiments in Astrophysics and Particle and Nu- 
clear Physics. Such data include: secondary par- 
ticle and isotopic production cross sections, total 
interaction nuclear cross sections and lifetimes of 
radioactive species, gas mass calibrations and gas 
distributions in the Galaxy (H2, H I, H II), inter- 
stellar radiation field (ISRF), CR source distribu- 
tion and particle spectra at the sources, and the 
Galactic magnetic field. All interactions that par- 
ticles might undergo during transport, such as en- 
ergy losses, and 7-ray and synchrotron production 
mechanisms, are similarly included. 

Study of the DGE will advance greatly with the 
forthcoming GLAST mission. In the following, we 
describe our ongoing efforts for understanding and 
modelling the DGE that will be incorporated into 
the model for the GLAST Large Area Telescope 
(LAT) Science Groups. 
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CR Propagation and GALPROP 

GALPROP is a code for CR propagation and dif- 
fuse 7-ray emission. We give a brief summary 
of GALPROP; for details we refer to the relevant 
papers [11, 14, 20, 21, 23] and a dedicated web- 
site^ . The propagation equation is solved numer- 
ically on a spatial grid, either in 2D with cyUn- 
drical symmetry in the Galaxy or in full 3D. The 
boundaries of the model in radius and height, and 
the grid spacing, are user-definable. Parameters 
for all processes in the propagation equation can 
be specified. The distribution of CR sources can 
be freely chosen, typically to represent supernova 
remnants. Source spectral shape and isotopic com- 
position (relative to protons) are input parameters. 
Cross-sections are based on extensive compilations 
and parameterisations |10|. The numerical solu- 
tion is evolved forward in time until a steady-state 
is reached; a time-dependent solution is also an op- 
tion. Starting with the heaviest primary nucleus 
considered (e.g., ^''Ni) the propagation solution is 
used to compute the source term for its spalla- 
tion products, which are then propagated in turn, 
and so on down to protons, secondary electrons 
and positrons, and antiprotons. In this way secon- 
daries, tertiaries, etc., are included. Primary elec- 
trons are treated separately. Gamma rays and syn- 
chrotron emission are computed using interstellar 
gas data (for pion-decay and bremsstrahlung) and 
the ISRF model (for IC). 

We are continuously improving the GALPROP 
code to keep up with new theory and data. Recent 
extensions to GALPROP relevant to the GLAST- 
LAT diffuse emission model include 

• interstellar gas distributions based on current H I 
and CO surveys (see below) 

• H2 mass calibration (Xco-factors) which can 
vary with Galactocentric distance 

• new detailed calculation of the ISRF (see below) 

• proper implementation of the anisotropic IC scat- 
tering [12] using the new ISRF (Figure 1 [left]) 

• new parameterisation of the tt" production in pp- 
collisions [9] which includes diffractive dissocia- 
tion 

• the extension of the 7-ray calculations from keV 
to tens of TeV, and the production of full sky maps 



as a function of energy; the output is in FITS for- 
mat (Figure 1 [right]) 

• a dark matter package to allow for propagation of 
the WIMP annihilation products and calculation of 
the corresponding synchrotron and 7-ray skymaps 

Interstellar Gas 

The maps of the neutral interstellar medium (ISM) 
used in the 7-ray intensity calculations have been 
updated recently. The neutral gas is traced by ob- 
servations of the 2 1 -cm line of H I and the 115 GHz 
line of CO (the standard surrogate for H2, which 
is not directly detectable at interstellar conditions). 
The differential rotation of the Milky Way causes 
distance-dependent Doppler shifts of the line fre- 
quencies. These shifts can be used to derive ap- 
proximate Galactocentric distances for the emit- 
ting regions corresponding to the observed spec- 
tral lines. We use the rotation curve of Clemens 
[3] in deriving Galactocentric distances and divide 
the Milky Way into equidistant rings of ~2 kpc 
width. Because for longitude ranges within 10° of 
the Galactic centre and anticentre so-caUed kine- 
matic distances cannot be determined, we interpo- 
late the maps across these ranges, using a method 
that ensures that the integrated column densities in 
the interpolated regions are consistent with the sur- 
vey observations. 

The new LAB survey of H I [8] is now used for cal- 
culating the 'rings' of H I. This survey has uniform 
coverage of the entire sky and has been carefully 
corrected for the effects of stray radiation. The 
CfA composite CO survey [4] is now used for cal- 
culating the CO rings. The data for the ring that 
contains the solar circle are augmented using a new 
intermediate and high-latitude survey; the initial 
results of this survey have been published [5], but 
additional observations have been made [6]. 

Interstellar Radiation Field 

The large-scale ISRF of the Galaxy is the result 
of stellar emission and dust reprocessing of the 
star light in the ISM. There is also a contribution 
by the cosmic microwave background (CMB). A 

L http://galprop.stanford.edu 
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Figure 1: Left: The ratio of anisotropic IC to isotropic IC for Galactic longitudes I = 0° and 180° vs. 
Galactic latitude [16]. Right: 7-ray spectrum of inner Galaxy (330° < I < 30°, |&| < 5°) for the optimised 
model. Vertical bars: COMPTEL and EGRET data, heavy solid line: total calculated flux. This is an update 
of the spectrum shown in [23]. 



model has been constructed for the Milky Way 
ISRF [15, 17, 18] incorporating details of stellar 
population distributions based on recent data from 
surveys such as 2MASS and SDSS, and a radia- 
tive transfer treatment of dust scattering, absorp- 
tion, and re-emission of the star light in the in- 
frared. The dust distribution in the model follows 
the gas distribution; to ensure the ISRF is consis- 
tent with the GALPROP code we use the gas dis- 
tributions described above. 

The ISRF model allows the calculation of the spec- 
tral energy density (SED) and angular distribution 
as a function of position and wavelength through- 
out the Galaxy. As an example of the model out- 
put, we show in Figure 2 the local SED (left) and 
the local intensity distribution at 2.2 fim (right). 
The SED is important for the CR electron energy 
losses during propagation. The intensity distribu- 
tion of the ISRF, which was previously not avail- 
able in the literature, allows the calculation of the 
IC emission using the anisotropic IC cross section 
[12]. This has been shown to produce significant 
differences over the sky when compared to the as- 
sumption of an isotropic ISRF (Figure 1 [left]); the 
latter approximation is true only for the CMB. 



Summary 

From the EGRET era we have learned a great deal 
about the DGE in the MeV-GeV range (see, e.g., 
[13] and references therein) while the GeV-TeV 
range remains largely unexplored. Recent VHE 
observations [1, 2] indicate the Galaxy is full of 
surprises. The DGE is also present at multi-TeV 
energies, but large variations are to be expected 
because of the inhomogeneity of the sources, and 
hence CR distribution; this contrasts with the case 
in the MeV-GeV range where the DGE has a sig- 
nificantly smoother distribution. 

We have given a brief summary of work that is be- 
ing done on the GLAST-LAT DGE model prior to 
launch; naturally, adjustments will be required af- 
ter launch. The GLAST-LAT wiU study the DGE 
in the GeV-TeV range, providing a clearer picture 
and connection between the spacecraft-borne and 
ground-based instruments. This will provide much 
new information on CR propagation and sources, 
and the ISM. 
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Figure 2: Left: Local ISRF spectral energy density. Line-styles: solid, total; thick dashed, stellar; thin 
dashed, scattered; chain, dust; dotted, CMB. Data points are summarised in [18]. Right: Model local ISRF 
intensity at 2.2 /im as a function of Galactic longitude I and latitude h. 
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